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Abstract 

We consider in a simple and general way elastic waves in isotropic and anisotropic 
media, their polarization, speeds, reflection from interfaces with mode conversion, and 
surface waves. Reflection of quasi transverse waves in anisotropic media from a free 
surface is shown to be characterized by three critical angles. 



I. INTRODUCTION 

In our time of supercollider, quantum computing, teleportation, dark matter and an 
eager search for a new physics, acoustics and elastic waves look like an archaic science. 
Am.J.Phys. rarely publish paper on this topic. In fact, we have found the single article 
of 1980 [1] relevant to our consideration. Indeed this science looks archaic, because every- 
thing seems to be well resolved and the science became like an engineering tool, frequently 
used in many applications. We will show that this impression is wrong. 

The theory seems to be well established (see, for instance [2]), and all the textbooks [3], 
are unanimous in its presentation. The main notion is a displacement vector u(r,t) of 
a material point at a position r at a time moment t. Its Descartes components Ui(r,t) 
obey the Newtonian equation of motion 

d 2 d 

where p is the material density, components of the radius-vector r, 0~ij IS cL stress 

tensor 

a ij = C ijklUkl, (2) 

which is proportional to the deformation tensor Uij 

1 / duj du k \ 
^ = 2(^ + ^-J' (3) 

and coefficients of proportionality in (2) comprise themselves a tensor with the sym- 
metries: 

C-ijkl C-jikl Cijlk Cklij- (4) 
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As is usual, in (1), (2) and everywhere below a summation over repeated indices is as- 
sumed. 

In isotropic media the tensor Cijki is very simple: 

Cijki = ^ij5 k i + fi(Su5 k j + S ik 5ij), (5) 

where Sij is the Kronecker symbol equal to unity for i = j and to zero otherwise, and A, 
ji are two parameters called Lame elastic constants. 

In the case of anisotropic media (usually crystals) the tensor contains considerably 
more parameters [4-7]. They are phenomenological, and their physical meaning is not 
sufficiently clear. 

The displacement vector u(r, t) is usually represented as a sum u(r, t) = Vy? + V x 1/? 
of two parts, where ip is a scalar, if) is a vector potentials, V is the differential vector 

d d d 

V = e x V x + e y V y + e z V z = e x — + e y — + e z — , (6) 

and e x ,y,z are unit vectors along three axes of an orthogonal Descartes reference frame. 

We will show that the use of the scalar and vector potentials is not necessary. They only 
complicate the theory. Everything can be presented in much more simple and transparent 
form with a wave function, like in particle physics. In the case of isotropic media such a 
presentation makes the theory looking almost trivial. 

In the case of anisotropic media it is quite instructive to consider not crystals, but a 
medium with an anisotropy distinguished only by a single direction [8], because a single 
vector is sufficient to elucidate the difference of isotropic and anisotropic media. We found 
very interesting counter intuitive features of waves reflected from interfaces in anisotropic 
media. Surprisingly, but these features have not yet been discussed in literature. 

After general introduction in the second section to the theory of elastic waves in 
isotropic and anisotropic media, we in the third section consider isotropic ones. In 
isotropic media elastic waves are very naturally break up into three classes, which are 
called modes: two modes with transverse and one mode with longitudinal polarizations. 

Theory of waves in a homogeneous medium is rather primitive. It becomes more rich 
when the medium contains an interface or a free surface. In that case a wave reflects, 
refracts and multiply splits at the interface. The presence of a surface or an interface gives 
rise also to waves of the fourth mode — the surface waves, which run along the surface, 
exponentially decay away from it, and have a mixed polarization. We consider how do 
they appear and what are their properties. 
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In the fourth section we go to anisotropic media with a single anisotropy vector. There 
again all waves break up into three classes-modes, but only one mode has purely transverse 
polarization. The other two are hybrids which are nor transverse nor longitudinal. One 
of these modes is called quasi transverse and the other one — quasi longitudinal, because 
the smaller is the anisotropy parameter, the closer is their polarization to pure transverse 
and longitudinal directions respectively. 

This is more or less evident and simple. The complications start when a wave meets 
an interface or a free surface. In that case we have reflection and refraction, which are in 
general accompanied with triple waves splitting. Specular reflection is absent, and surface 
waves become exotic. An experience acquired from isotropic media considerations leads 
to a conclusion that in anisotropic ones an incident wave can completely transform into a 
surface wave! It is absolutely unacceptable because of the energy conservation law. The 
incident wave carries an energy, which must accumulate in the surface one, therefore the 
amplitude of the last one must grow exponentially. It is impossible to describe such a 
process by a linear stationary wave theory. So we meet a paradox: the theory predicts 
some non physical solutions, which cannot be described by the theory. 

We considered the problem in details, and found that the reflected waves do not ac- 
cumulate into a surface mode because of unexpected counter intuitive properties of the 
elastic waves in anisotropic media. For instance, besides the first critical grazing angle tp c , 
similar to that one in isotropic media, at which a quasi longitudinal wave becomes of the 
surface type, there is a second critical angle, (p cl < ip c , when the reflected wave, which is 
naturally thought of as moving away from the reflecting surface, changes its direction, as 
if starting to move toward it, though its energy flux remains going away from the surface. 

There is also a third critical angle, y? C 2 < ip c i, at which the energy flux toward the 
surface becomes zero. This angle is not zero, i.e. <p C 2 > 0, and annulation of energy 
flux at this angle means that we cannot direct a ray of elastic waves to the surface at an 
angle (p < (p c2 . It is really strange, but we consider it as a hint, which the linear theory of 
elasticity gives us to point out the cases, where we have to involve a nonlinearity. We could 
not understand how to introduce it, but we were so much impressed by the unexpected 
properties of the elastic waves, that we decided to relate about them to the readers of 
this journal. 

Our research was started because of a need to explain the difference between theo- 
retically predicted and experimentally measured anisotropy of sound speed in rocks [8]. 
The theory, based on texture of rocks measured at a neutron diffractometer, predicted 
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anisotropy of speeds in some rocks three times lower than the one measured with ultra- 
sound. We hope that our results will shed light on this difference, and will be checked in 
an experiment of the kind discussed in [1]. 



II. THE MAIN EQUATIONS FOR ELASTIC WAVES 

The starting point for study of elasticity is the free energy density of a medium defor- 
mation [2]. For isotropic media it is 

F = ^ul + nul, (7) 

where tty is the deformation tensor (3), A, \i are the Lame elastic constants, notation ufj 
means uijUji, and summation over repeated indices is assumed. 

In anisotropic case we have to distinguish a direction, say along a unit vector a, and 
introduce a new elastic constant, say (, which is of the same dimension of energy density 
as A and /i. Then the free energy becomes 

F = + ~ ([(ajUji) 2 + (ujiai) 2 }, (8) 

where a,j are Descart components of the vector a, anisotropic part is proportional to 
square of the vector with components djiiji, and the sign before ( is not necessary negative. 
Though the two anisotropic terms in square brackets are identical {u\j is symmetrical) we 
put them separately to get the tensor with required symmetry (4). 

The anisotropy term shows that the energy of deformation depends not only on change 
of volume (the term oc A) and change of shape (the term oc fj,), but also on angles of both 
deformations with respect to the anisotropy vector a (the term oc (). 

A question can be raised here: why do we use such an anisotropy modification of the 
isotropic free energy (7)? Is it not possible to find a different one? The reply is yes. It 
is possible to use a different modification. For instance, instead of (8) we can accept the 
free energy in the form 

F = ^ u u + H u % - C( W a ') 2 > ( 9 ) 
or we can use a combinations of anisotropic terms in (9) and in (8). We chose (8), and it 
is an arbitrariness. We think that consideration of different types of anisotropy is a good 
task for students. 

With the free energy we can define the stress tensor 

dF 

= < 10) 
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which after substitution of (8) gives 

Gij = XSijUn + 2/jiUij - 2((a i u jl ai + a^ H %) = c ijkl u kl , (11) 

where 

(Hjki = X5ij5 k i + n(5u5 k j + S ik 5ij) - C(ai5jia k + 5uaja k + a { bj k ai + 5 ik ajai). (12) 

The expression for isotropic medium is obtained in the limit ( — > 0. Note that the tensor 
Cij k i satisfies the symmetry requirements (4). In the case of (9) the anisotropic part of 
the tensor c^ki will be simply a product (ai(ija k ai. 

With the stress tensor (11) the Newtonian equation of motion (1) for the displacement 
vector becomes 

piii = y j(Tij = p[A Ui + Vi(V • «)] + AVi(V • «)- 

-C (oi[A(u ■ a) + (a • V)(V • «)] + (a • V) 2 ^ + V^a • V)(o • «)) , (13) 

where p is the medium density. 

We can seek solution of (13) in the form of a complex plain wave u(r, t) = A exp(ikr — 
iut), like a wave function in particle physics, where vector A ia a unit polarization vector. 
Of course, the elastic waves are real waves, so they are represented by the real part of 
the complex wave function. Later we shall see, where we must be especially careful in 
description of elastic waves with complex function, but for now we meet no difficulties, 
and after substitution of such u into (13) we obtain an equation for A: 

pu 2 A = pk 2 A + (A + p)k(k ■ A) — 

-C (a[k 2 (a ■ A) + (k ■ a)(k ■ A)] + (k ■ a)[(k ■ a) A + fe(o • A)]) . (14) 

It is convenient to transform this equation to dimensionless form dividing both parts of 
(14) by pk 2 . After introduction of the standard transverse speed c t = \J p/ p, the phase 
speed of the wave V = cu/k, dimensionless ratio v = V/c t , dimensionless parameters 
E — (A + p)/p, £ = C/p and the unit vector n = k/k the equation becomes 

tt 2 A = Ek(k ■ A) - f [{a(a ■ A) + (k ■ a) cl(k ■ A) + A(a ■ k) + k(A -a)]}, (15) 

where Q 2 = v 2 — 1. First we put £ = 0, and consider an isotropic medium. 
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III. WAVES IN ISOTROPIC MEDIA (f = 0) 

For isotropic media the equations (13-15) are reduced respectively to 

put = V jUiJ = fi[Au t + Vi(V • «)] + AV,(V • u), (16) 

pu 2 A = fik 2 A+{\ + n)k(k- A), (17) 

Vl 2 A = Ek(k ■ A). (18) 

For the given propagation direction k we can introduce two orthonormal vectors and 
e( 2 \ which are perpendicular to k. In the orthonormal basis e^, e^ 2 \ k the polarization 
unit vector A is representable as 

A = a^e^ + a^e^+p K . (19) 

After multiplication of (18) consecutively by e^ 1,2 ) and n we obtain three equations for 
coordinates a^ 1 ' 2 ) and (5: 

ft 2 a (1 ' 2) = 0, [VI 2 - E]f3 = 0. (20) 
They are independent and give three solutions 

A (l,2) = e (l,2) > A (3) = K (21) 

Their speeds are determined by equations 



( v (1.2))2 = 1 ^ V (l,2) = ^=Ct, (V^) 2 = E+l - = V^+lCi = ^/ = Q- 

(22) 

Since V = u/k, we can tell that for a given frequency u the wave numbers k — \k\ of 
three modes are 

/c (1 ' 2) = -, /c (3) = -. (23) 

The two unit vectors e^ 1 ' 2 ) above are orthogonal to each other, and lie in the plane, per- 
pendicular to the unit vector k, but their azimuthal angle around k can be arbitrary. We 
can use this freedom to facilitate solution of different problems. In particular, below, when 
we consider reflection from an interface. There we can choose to be perpendicular to 
the incidence plane, and to be inside it. 
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A. Reflection from an interface 



Suppose that the medium consists of two parts: one at z < with constants A, /i, p 
and another one at z > with constants A', //, p', then there is reflection and refraction of 
waves at the interface z — 0. If a plane wave u inc {r, t) = Aexp(ikr — iujt) incident from 
z < is of mode (j is one of the numbers 1, 2 or 3) then the interface transforms 
this displacement vector to 



u(r, t) = exp(ife|iT*|| — iujt) x 



x 



A w e *¥> + ^ r o-o A g e -^ 2 ) S(z < 0) + J2t {jl) A^e ik '" >z e(z > 0) 



■„(0 



'(0. 



i=i 



i=i 



(24) 



where r^ l \ ft 1 ) are reflection and refraction amplitudes of the l-ih mode (/ = 1,2,3) for 
the incident j-th mode, and 6 is a step function, which is equal to unity, when inequality 
in its argument is satisfied, and to zero otherwise. 

Below the interface (z < 0) the displacement consists of the incident wave of j-th 
mode and reflected waves of modes (the lower index R means reflected). Above the 
interface [z > 0) the displacement consists of transmitted waves of modes Aj) (the lower 
index T means transmitted). 

All the waves differ from one another not only by polarization, but also by the wave 
vector k, which can be represented as 



k = fey + k± = rkn ± nk±, 



(25) 



where r is a unit vector along the interface, and n is a unit vector along z-axis perpen- 
dicular to the interface, as is shown in Fig. 1. Note, that the component, k\\ = rk\\, of 
the wave vectors is identical for all the waves, as is demonstrated by the first common 
factor exp(ife||T*|| —iut) in (24), where r\\ are coordinates in the interface. The vector k\\ is 
identical, because the space along r is uniform and nothing can change this component. 

The normal components and k'jp of wave vectors of modes A® and Aj) respec- 
tively, are positive numbers and their value depends on the mode /. Since k\ = k 2 — k 2 , 
then from (23) it follows that 



— -k 2 



k 



'(1,2) 



— ~k 2 



k 



(3) 



^_ _ 1.2 



fc?> = 



\^2~ k l ( 26 ) 



where c t j, and c' tl , are the speeds defined in (22) for lower and upper spaces respectively. 
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To find reflection and refraction amplitudes, we need boundary conditions. One of 
them is continuity of the displacement vector: 

3 3 



«|z=-o — u\ z=+ o 



A^ + Y. r{3l)A R=Y. t(jl)A T, (27) 
i=i i=i 

and the second one is the continuity of the stress vector T with components Tj = cr^raj. 
According to (11) this vector for a displacement u(r,t) is equal to 



T(u(r, *)) = An(V • u) + n\V{u ■ n) + (n • V)u]. 



Continuity of the vector T is equivalent to the equation 

3 3 



(28) 



(29) 



where the vector £? is defined as 



B = -iexp(-ifcr)T(Aexp(ifcr)) = An(fc • A) + fi[k(A ■ n) + A(n ■ k)], (30) 

for every plane wave Aexp(ifcr). 

Note that the condition (29) makes it possible to continue the wave equation (13) from 
z < to z > 0. If it is not satisfied, the differentiation of in (13) creates 5(z)-function 
and the wave equation becomes inhomogeneous [9] with a source term at the interface 
z = 0. 

To find reflection and transmission amplitudes we need to multiply both equations (27) 
and (29) by three mutually orthogonal unit vectors to get in general 6 equations for 6 




FIG. 1: Reflection of a transverse wave from an interface between two different isotropic 
media. Reflected and refracted waves contain two modes: A^' 3 ^ and A^' 3 ^ respectively. The 
A^i mode goes at specular grazing angle (p^ = ip, the longitudinal mode A^ goes at grazing 
angle (p^ < <p( 2 \ 
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unknowns. It is convenient to choose the right triple the vectors r, n, as shown in 
Fig, 1, where the vector e® is perpendicular to the incidence plane and in Fig. 1 points 
toward the reader. 

Before calculations of the amplitudes of all the reflected and refracted waves we can 
easily understand what are all the angles. The grazing angle y? of the wave with the wave 
vector k is defined via relation cosy? = r ■ k/k = k\\/k = k\\V/uo. Since k\\ and u> are 
identical for all the waves therefore the value cos ipk>) /V^ are also the same for all the 
waves. And because of (22) we can write 

cosy? cosy?' 1 ' 2 ) cosy?® cosy?^ 1 ' 2 ) cosy?'® . . 

where y?W <p'W denote grazing angle for respectively reflected and transmitted waves of 
mode i, and y? without indices denotes the grazing angle of the incident wave. 

To find directions of polarization after reflection is very easy and we leave it as an 
exercise for the reader to check that with account of (23) 

/)_d2) xe (i),__^Il¥ ,0) -kfn + k\\r 

1. Reflection of A® mode 

The simplest is reflection and refraction of A® mode. Its polarization is e®. After 
multiplication of equations (27) and (29) by e® we get 

1 + r (n) = ^(1 _ r ^)k ± = fi 2 t^k' ± , (33) 

from which it immediately follows that 



,(n) 



fj,k± — /J,'k'j_ 
/ik± + fi'k' ± ' 



(34) 



where k± = ^Ju 2 /c 2 — A; 2 , and k' ± = ^Ju 2 /c' 2 — k 2 . We see that this mode is reflected 
specularly and no other modes are created. 

2. Reflection of A® mode, and the mode conversion. 

The more interesting is the case of the incident A® mode shown in Fig. 1. Its 
reflection and refraction creates longitudinal mode A®, and because the speed q of A® 
is larger than the speed c t of the specularly reflected A® mode, the grazing angle y?® is 
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less than p( 2 \ Therefore we can expect that at some angle ip — ip c of the incident wave, 
the angle (p^ becomes zero, which means that the longitudinal mode ceases to propagate 
in the direction z < 0. Since according to (31) cosp^ = (ci/c t ) cosp < 1, we find that 
ip c = axccos (cf/cj). The similar considerations are applicable to the refracted waves, and 
we can expect that at p < arccos(max(ct/q, c' t /c'j) there appears a purely longitudinal 
surface wave propagating along the interface. 

To find amplitudes of the reflected and refracted modes we have to multiply the two 
equations, (27) and (29), by n and r. As a result we get a linear system of four equations 
for 4 unknown A 22 \ A 23 \ and t@ 3 \ which can be solved analytically. However it is a 
boring job, so it is better to pass it to computer. 

The analytical solution can be found for reflection from a free surface, where we have 
a single boundary condition 

B (2) + r (22) B g) + r (23) B g) = Q (35) 

In this case we have only two reflected waves and multiplication of (35) by n and r with 
account of (30) and (32) gives only two equations 

Lll {l-r^)+r^—- 3=0, (36) 



W) " (! + r {22) ) - 2r (23) ^ = 0. (37) 



Their solution is 



r (23) _ ^ 4fcj 2) fe||(A:( 2 > 2 - 2fc^) _ 4fcj 3) fcj 2) fc} - (k^ 2 - 2fcj) 2 

Mfhfh 2 + (k^ 2 - 2k\) 2 ' ~ Akfkfh 2 + (A;( 2 ) 2 - 2A;J)2' (3§) 

At </? = (p c the vector k^ 3 ) of the longitudinal wave propagation direction coincides with 
r, and therefore the length k^ = u/ci of the wave vector becomes equal to k\\. When 
the grazing angle (p decreases below p c the value of k\\ increases, but tv/ci does not change. 
Therefore at p < p c we get k\\ > uj/q, and k^ = ^Juj 2 /c 2 — k 2 becomes imaginary. We can 
denote it —iKi. With such a normal component of the wave vector the longitudinal wave, 
propagating along the free surface, becomes localized in the layer of thickness I = l/Ki, 
where K x = ^jk 2 — uj 2 /cf. In other words, it becomes longitudinal surface wave (lower 
index S means surface) with complex polarization vector A s . But what a strange wave 
it is! Since the incident wave can have arbitrary u and ku < k — uj/c t , the longitudinal 
surface wave with the same u> and k\\, has the speed along the surface, (denote it Vg ) 
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equal to V<j — uo/k\\ > c t , which means, that it is not the Rayleigh surface wave, because 
the speed cr of the Rayleigh wave, as is well known, is less than c t \ 

However, really, it is not strange. This longitudinal surface wave satisfies the same 
wave equation Q 2 = E of (20), and has longitudinal polarization 



4 3) = (fc|| T - iKmyyjq + Kf, 



(39) 



i.e. its normal component is imaginary. It is not dangerous that this polarization is a 
complex vector. The displacement must have a real value, therefore the displacement 
with a complex polarization vector is 



u^p oc Re 



(k\\T — iKin)e 



ikuru+Kiz— iwt 



(fc||TCOs(fe||r|| — uot) + Ifjn sin(fc||r|| — uot) 



e K ' z , 
(40) 

i.e. the phase of oscillations along vector n is shifted by n/2 with respect to oscillations 



(3) 

along vector r. The speed of the longitudinal wave along the interface, V s — u 
be arbitrary, though because of fcjj > uj 2 /c 2 , this speed lies in the interval c t < Vg 3 ^ < q. 

In a similar way we can define the transverse surface wave. It satisfies the wave equation 
Q 2 = of (20), and has transverse complex polarization 



Ail. 



can 



A (2) 



l{ n + iK t r)/Jk 2 + K 2 , 



where K t = Jk 2 — u 2 /c 2 . The real displacement vector in it is 



u§ ^ oc Re 



(k\\n + iK t r)e 



ikara+Ktz—iut 



(/c||ncos(fc||7 , || — uot) — _K" t rsin(fc||r , || — uot) 



(41) 



e KtZ , 
(42) 



(2) 



(2) 

and the speed along the interface can be arbitrary but less than c t . 

When longitudinal wave is of the surface type, the reflection amplitude of the 
mode according to (38) becomes 

(k^ 2 - 2kf,) 2 - AiKikfhf, 



r (22) = 



(43) 



{k^ 2 -2k 2 ) 2 + AiK l k { l ) k\ 
It is a unit complex number, therefore it describes the total reflection of the incident wave. 



3. Energy flux distribution between two reflected waves 

Because of energy conservation the energy flux density of the incident wave along the 
normal to the interface must be equal to the sum of energy flux densities of the reflected 
waves. Let's check, whether they are really equal. 
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When the displacement and therefore the stress tensor are real functions, the energy 
flux density of an elastic wave is described by a vector j with components 



ji = -(a u dui/dt), 



(44) 



where (F) means averaging of the function F over time. We use displacement in the 
form of complex plane waves, but the energy flux density should have only a real value, 
therefore the Eq. (44) can be represented as [10] 

: dui 



Ji 



*dm du^ 



= -Re 



dt 



(45) 



where Re(F) means real part of F, and * means complex conjugation. We are interested 
in the flux density along the normal n to the interface, therefore we need to calculate 
nj = ujR,e(iT*u), where we used (28). Taking into account the definition (30) we can 
represent the energy flux in the form 



= -Re{B* ■ A) = ERe((A* ■ n){k ■ A) 



k 



In the case of the incident A^ mode its incident flux is = pwk 



(2) 



(46) 



puc 2 t k ( l ] . 



Reflected fluxes of the two modes are j$ = pcflr^l 2 ^ , j$ = pcf |r^ 3 ^| 2 A;^ 3 ' ) Energy 
conservation law requires 



cfk^ 



I (22) 1 2 , I (23)|2 W "a _ -, 



(47) 



Substitution of (38) shows that this equation is satisfied. 

When the longitudinal wave becomes of surface type, it does not produce a flux from 
the surface. Therefore the law of energy conservation (47) reduces to 



,(22) I 2 



1, 



(48) 



which, according to (43), is also satisfied. 



B. Energy density of the longitudinal surface wave 

There is also one interesting question: what the energy density is accumulated in the 
longitudinal surface wave. This question is interesting, because it is this wave can be 
important for predictions and estimation of magnitudes of the earthquakes. 

The vector of energy density flux of a wave of mode A according to (45) can be 
represented as 

J = Re(pujc 2 t ul[EA*(k* ■ A) + k*]), (49) 
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where we introduced an amplitude uq of the wave. Expression (49) is valid for real and 
complex wave vectors and polarizations. The absolute value of this flux for the incident 
transverse wave is 

J = pucfulk. (50) 
For the surface longitudinal wave with account of its factor r 23 we get 

jf (z) = pujul\r 2 *\ 2 c 2 k\\ exp(2iT,z). (51) 

We do not know u , so we can find only ratio Q = Jg(z = 0)/J. Substitution of r 23 from 

10 
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FIG. 2: Dependence of Q on cos ip calculated for E = 1.2. The left point on the abscise axis 
corresponds to critical cos(p c = 0.674, where Q = 323.16. 



(38) gives 



Q 



k- - 2 



or 



fc(2) 



Q = 



SL\r 23 \ 2 



16fci 2) 2 k 2 (k^ 2 - 2k 



2\2 



lQK 2 k { l )2 kl + (k^ 2 - 2k\Y'' 
cosy? sin 2 (4<p) 



(52) 



(53) 



4(cos 2 <f — cos 2 <f c ) sin 2 (2(/?) + cos 4 (2(/?) 
Dependence of this function on cosy? is shown in Fig. 2. We see that the highest energy 
density is accumulated in longitudinal surface wave, when ip is slightly less than ip c . There 
is also a maximum at small angles where the ratio Q is close to unity. 



C. The Surface Rayleigh wave 

We considered above the two surface waves, which satisfy the wave equations, but 
cannot exist independently, because without the incident and reflected waves they do not 
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satisfy the boundary condition. The Rayleigh surface wave exists without the incident 
one, and its speed cr — uj/k\\ < c t is fixed. To get equation which determines this speed 
Cr we represent the boundary condition (35) in the form 

^« + B<? + ^Bg»=0, (54) 

where r^ 22 ) and r^ 23 ) are given by (38). With (54) we can immediately find the speed of 
the Rayleigh surface wave. It corresponds to such ou/k\\, for which the first term in (54) 
is zero. Since ^ 0, therefore the first term is zero only when 

1 



,(22) 



= 0. (55) 



In such a case the incident wave disappears, and the whole wave field contains only two 
waves propagating along the free surface. 

Let's remind that a similar trick helps to find bound states of particles in quantum 
mechanics. Reflection of a particle from a one dimensional potential well is described in 
asymptotic region x — > — oo by the wave function exp(ikx) +r(k) exp(-ikx), where r{k) is 
a reflection amplitude. This wave function can be also represented as (l/r(/c)) exp(ikx) + 
exp(— ikx). In bound states the wave function at x — > — oo has asymptotics exp(— Kx), 
where — K 2 is proportional to the bound level E^. To find K we need to solve equation 
l/r(k) = 0, which annuls the incident wave exp(ikx). Every root of this equation k n = 
—iK n corresponds to n-th bound level E bn oc —K 2 . In that respect the Rayleigh surface 
wave is a bound state of elastic waves, 

After this digression we go back. From (43) it follows that (55) is satisfied, if 

4/fci 3) A:i 2) A:f + (A:^ 2 > 2 - 2/fcf) 2 = 0. (56) 

It is important to note that the third term in (54) does not disappear though it also 
contains the factor 1/r^ 22 ^. It does not disappear because and according to (38) 
have the similar denominators, and they cancel each other in the ratio r^ 23 ) / >( 22 ). In fact 
the amplitudes and r^ 22 ) play equal roles, so instead of (54) we can write 

1 ^(22) 

and seek solution of the equation l/r^ 23 ) = 0. The result will be the same. 

Let's denote the speed of the wave propagation, w/k\\, along the interface by cr (speed 
of the Rayleigh wave), and its ratio to c t by x = Cft/c t . Since k]_ and k ± in surface 
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waves are to be imaginary then k]_ = uj 2 /c 2 — k 2 < 0, k]_ = uj 2 /cf — k 2 < 0, and the 
equation (56) is reduced to 

4v / T^ : Vl - ^x 2 = (2 - x 2 ) 2 , (58) 

where q = c t /ci. This equation and its solution are well known and can be found in all 
the textbooks on elastic waves. It can be solved even analytically for arbitrary q, because 
it is equivalent to an algebraic equation of 4-th order with respect to variable z — 1 — x 2 . 
When (58) is satisfied, the total displacement vector of the Rayleigh wave becomes 

(2) , r(23) (3) ft-ns 

U R (XU S + ^2j u s , (59) 
Substitution of (38), (40) and (42) into (59) with the same oo/k\\ = cr gives [11] 

u R (xncos(k\\r\\-ujt)[2q tqi e KlZ - (1 + q 2 t )e KtZ ] - r sin(Jfe||T- N - Lut)q t [2e K ' z - (1 + q 2 t )e KtZ }. 

(60) 

where K u = k^l-c 2 R /c 2 tp and q t>l = K u /k\\ = ^l-c 2 R /(? tl . 

In a similar way we can find the Stoneley surface wave propagating along the interface 
between two isotropic media. Though there are no principal difficulties, we do not consider 
it here because of technical complications. 

From (30) we can immediately conclude that the surface waves with polarization along 
the surface and perpendicular to direction of propagation do not exist, because the conti- 
nuity of the stress vector B requires continuity of the normal derivative of the displacement 
vector, which cannot be satisfied. 



IV. WAVES IN ANISOTROPIC MEDIA (£ / 0) 

In isotropic media it was natural to describe polarizations in an orthogonal basis e^, 
e^ 2 ) and k — k/k, which constitutes the right hand triple of unit vectors. The choice of 
and had some freedom because these two vectors can be rotated by an arbitrary 
angle around n. It was only in studying of reflection from an interface, where orientation 
of was fixed by the plane of incidence. In anisotropic media besides k we have 

also vector a, so for orientation of vectors e^, it is better to take the plane of vectors 
k and a into account, choosing in the plane, and = [e^ x k] perpendicular to it. 

After substitution of (19) into (15) and multiplication by three vectors and 
k = k/k we obtain the system of linear equations for c^ 1 ' 2 ), and (3. It has a solution, if its 



15 



determinant is equal to zero. This condition gives three possible possible speeds V^(k) 
(i=l,2,3) for the three wave modes. 

The simplest equation is obtained after multiplication of the (15) by The result 

is 

fiV^+^-ajV 1 ' = 0. (61) 

This Eq. is equivalent to 

VI 2 a) 2 = 0, (62) 

and it gives the speed of this transverse mode 

V {1) = u/k = -£(/«• a) 2 = ctyfl -£cos 2 #, (63) 

where 9 is the angle between vectors k and a. We see that this speed is less than c t , and 
it changes with 9. 




FIG. 3: The transverse, quasi transverse, VV\ and quasi longitudinal, V^ 3 \ speeds of 

elastic waves in an anisotropic media in dependence on cos#, where 9 is the angle between 
wave vector k and anisotropy vector a, for E = 1.5 and two different anisotropy parameters a) 
£ = 0.5; b) £ = 0.7. The unity on ordinate axis corresponds to q = yV/P- 

After multiplication of (15) by and n we obtain a system of two equations 

(n 2 + 0« (2) + 2£(a • e^)(K ■ a)P = 0, (64) 

[VI 2 - E + 4£(k • a) 2 ]f3 + 2£(a • e (2) )(^ • a)« (2) = 0, (65) 

where in (64) we used relation (a-e^) 2 + (n-a) 2 = I. We see that polarizations along 
and k are not independent. They combine and create two new hybridized polarizations, 
which we call quasi transverse and quasi longitudinal modes and denote by A*- 2 ' 3 -* like in 
isotropic case. 
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The system (64,65) has solutions, if 

[n 2 + f ] [n 2 - E + 4£(k • a) 2 } - A£ 2 (k ■ a) 2 (a ■ e^) 2 = 0, (66) 
From which it follows 

2(ft (2 ' 3)2 +S)=E + £[1- 4(k • a) 2 } t V '{ E + Cl 1 " 4 ( K • «) 2 ]} 2 + 16f 2 (« ' «) 2 (« ' e(2) ) 2 - 

(67) 

Since (k • a) 2 = cos 2 0, and (a ■ e^ 2 )) 2 = sin 2 9 then, because fi( 2 ' 3 ) 2 = \/( 2 ' 3 ) 2 /c 2 — 1, we 
get that the equation (67) is equivalent to 



= c t 

At £ — > their values are 



£ + £(1 - 4 cos 2 #) =F y [-E + £(1 — 4 cos 2 0)] 2 + 4£ 2 sin 2 (20) 
\ 1 £ + 2 • ( 68 ) 



V (2) « ct(l - - 0(£ 2 ), V< 3 > » c, - 2| cos 2 9 + 0(e), 



(69) 



where 0(£ 2 ) denotes a small number proportional to £ 2 . So can be called quasi 
transverse and — quasi longitudinal speed. 

All the speeds V^ l \ and depend on angle 9. This dependence is shown in 
Fig. 3. We see that if the anisotropy parameter £ is sufficiently large some modes at small 
angles 9 cease to propagate, because their speed, as is shown in Fig. 3b) for quasi transverse 
mode, does not exist. This speed becomes imaginary, therefore the wave number of the 
mode, = lu/V^ 2 \ also becomes imaginary, and the wave does not propagate. Of course 
it corresponds to too large anisotropy parameter. Since E — 1.5 then £ = 0.7 means that 
anisotropy energy ( is larger than the Lame index A, and in some directions the higher is 
deformation the less is the elastic energy, which is nonphysical. For smaller £ the speed 
V 2 is at no angle imaginary. 

From (64) and (65) it follows that polarization of propagating quasi transverse, 
and quasi longitudinal, modes are 

A(2) = £sin(20)e( 2 )-(^ 2 ) 2 +Q^ A(3) = (^ 3 ) 2 + Q/,-£sin(2g)e( 2 ) 
^(fi^ + O' + e 2 sin 2 (20)' v /(fi( 3 ) 2 +O 2 + e 2 sin 2 (20)' 

At small £ they, as can be expected, are: 

A (2) w e ( 2) + I sin ( 2 #) K) A (3) w K _ I sin ( 2 #) e ( 2 ). (71) 
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FIG. 4: Splitting of reflected and refracted waves at an interface between two different anisotropic 
media when the incident ia a wave of quasi transverse mode A^ 2 \ S denotes the specular 

(2) 

direction. Anisotropy vector a has such a direction that the speed of reflected A R mode is 
higher than that of the incident one. 

A. Reflection from an interface 

Reflection of waves from an interface in anisotropic media is in general characterized 
by trirefringency, as was correctly pointed out in [1]. An incident wave at an interface in 
general splits into three reflected and three refracted waves, and no wave is reflected spec- 
ularly. In Fig. 4 we present the scheme of reflection and refraction of a quasi transverse 
wave from an interface between two anisotropic media with different anisotropy vectors 
a and a' and different parameter p, A, /i and (. The anisotropy vectors in general are 
not in the incidence plane. In Fig. 4 they are inclined down, so the reader sees their tails 
denoted by crosses. 

The grazing angles <p®, <p' {i) (an gles between wave vectors k RT of reflected and re- 
fracted modes A^ T and the unit vector r) in the case when the incident wave is of mode 
j, are determined from the relations equivalent to (31): 

cos ip cosyW cosyj'w 
v R v T 

The value of the speed of a wave depends on the angle 9 between the direction of 
its propagation and the anisotropy vector a. It may happen that after reflection all the 
speeds are higher that the speed of the incident wave. Then the grazing angles of 
all the waves become less that that of the incident one as is shown in Fig. 5a), and we 
can expect that at some critical angle ip = ip c all the reflected and transmitted waves will 



18 




FIG. 5: Reflection, refraction and splitting of waves at an interface between two different 
anisotropic media, when the speeds of refracted waves are higher than that of the incident 
quasi transverse wave, a) The grazing angle of the incident wave is sufficiently high, and all the 
created waves are able to propagate in z direction, b) Unreal situation, when the incident quasi 
transverse wave is transformed into the surface wave containing all the three modes. 

accumulate into a single surface wave as is shown in Fig. 5b). 

Physically such a result is unacceptable, because the incident plane wave gives the 
energy flux toward the interface, therefore the energy must accumulate in the surface 
wave and the surface wave amplitude should increase with the time exponentially. We 
are dealing with stationary waves, therefore exponentially growing functions are excluded 
from our solutions. 

We should look what is wrong in our logic, considering an example, in which everything 
can be solved analytically. The analytical solution can be found in the case of reflection 
of a quasi transverse wave from a free surface, when anisotropy vector lies in the incidence 
plane, as is shown in Fig. 6. In this case we have only two reflected modes: quasi transverse 
and quasi longitudinal ones, and to find their reflection amplitudes we need to solve only 
system of two linear equations. 

B. Reflection of quasi transverse wave from a free surface, when anisotropy 

vector is in the incidence plane 

Let's consider reflection of a plane wave of quasi transverse mode from a free 
surface, when the anisotropy vector has such a direction, that the reflected quasi transverse 



19 




FIG. 6: Reflection of a quasi transverse wave from a free surface, when anisotropy vector lies 
in the incidence plane, and the reflected wave speed is higher than that of the incident one. 
Reflection is accompanied with creation of quasi longitudinal wave, a) The grazing angle of the 
incident wave is sufficiently large and both reflected waves can propagate in z direction, b) The 
grazing angle of the incident wave is sufficiently small, and quasi longitudinal wave propagates 
only along the surface. 

wave has higher speed than the incident one. 

The angles of reflected waves are determined by (31) 

cosy? cosy^ 2 - 1 cosy/ 3 ) 

where 6 and #( 2,3 ) are the angles between a and directions of propagation k of the incident 
and k^ 2 ' 3 -* of the reflected waves respectively. In these equations we do not know (0®), 
therefore we cannot directly find (p^\ Instead we have to use these equations to find both 
</?W and V ir {0^) simultaneously. 

Let's denote a = Tcosy? a + nsiny? a , n = rcosy? + rising and k' 1 ' = Tcos^p^ — 
nsinyW then cos# = a ■ k and cos0^' = a ■ Substitution of cos# and cos 9^ into 
(68) and after that into (73) gives transcendent equations that can be solved numerically. 

The result of calculations for cosy? a = 0.4, £ = 0.4 and E — 1.5 are shown in Fig-s 7 
and 8. Fig. 7 shows dependence of all the speeds on cosy?, and Fig. 8 shows dependence 
of cosy?'- 4 -' on cosy?. 

We see that at cosy? > 0.92 no reflected wave can propagate. What does happen there 
is the most interesting question! 

At cosy? < 0.5, where both reflected waves do really exist, we can find their reflection 
amplitudes. For that we have to solve the boundary condition equation 

B (2) + r (22) B (2) + r (23) B (3) = , (74) 
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COS((|)) 



FIG. 7: Dependence of speeds of the quasi transverse incident, V^ 2 \ quasi transverse reflected, 

(2) (3) 

V£ ' and quasi longitudinal reflected speeds on cos y? of the grazing incidence angle, when 
(a • t) = 0.4; £ = 0.4 and E = 1.5. We see that in some range of cos y> both reflected speeds are 
higher than that of the incident one. 

where B is defined like in (30): 

B = -iexp(-ikr)T(^Aexp(ikr)^j = (n ■ k)A + k{n ■ A) + (E - l)n(A ■ k)- 

-£{a[{n ■ k){a ■ A) + (a ■ k)(n ■ A)] + (n ■ a)[{a ■ k)A + k(a ■ A)]}. (75) 

Multiplying (74) by n and r, we obtain 2 equations for two reflection amplitudes. Let's 
denote p = (n • B (2) ), 5 = (t • B (2) ), & = (n • B§ } ) and <J f = (r • B$) where i=2,3, then 
Eq. (74) can be represented in matrix form 

cm: a q-* 

and its solution is elementary. We do not represent the final analytical result because it 
does not look sufficiently informative. 

u 



CQS(0< 2 >) 

( 

°0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

oos(<|>) 

FIG. 8: Dependence of cosy/ 2 ) and cosy/ 3 ) on cos <p. We see that both cosy/ 2 ) and cosy/ 3 ) 
reach unity at cosy? < 1. It looks as if both waves become of surface type, when the incident 
one still remains to be the plain wave. 
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FIG. 9: Dependence of reflectivities R( 22 ^ = 1 7- ( 22 ) | 2 of quasi transverse and = |r( 23 )| 2 of 

quasi longitudinal waves on cos f>. b) The difference of energies AE of the incident and reflected 
waves. We see that all the values are well calculated only up to cos ip = 0.5, where longitudinal 
wave becomes of the surface type. 

The numerically calculated dependence of reflectivities = |r^ 22 ^| 2 and = 

I r (23) 1 2 on COSL p j s shown in the left panel of Fig. 9. The correctness of calculations is 
supported by the panel b), which demonstrates the energy conservation at reflection. All 
the calculations are possible only up to cosp = 0.5. Above this value the the wave vector 
of the quasi longitudinal wave becomes complex, and equation (73) has no solutions. 



C. Waves propagation near a free surface 

To see what happens, above cosip = 0.5 we have to make calculations differently. At 
the interface there are two conserved values: the frequency to and the wave number kn 
along the interface. It is worth to divide both parts of the Eq. (14) by fik 2 , to introduce 
the value 

„ ^ V^ 2 (cos6) 

T = = = -, (77) 

c 2 k 2 k 2 c 2 cos 2 <f c 2 cos 2 f 

and the dimensionless wave vector k = k/k\\ = t + qn, where q = k±/k\\. After that the 
Eq. (14) is transformed to 



T-l-q 2 + ^k-a) 



A = Ek(k ■ A) 



-£ (o[(l + q 2 )(a -A) + (k- a)(k ■ A)] + k(k ■ a) (a ■ A)]) . (78) 

This equation describes propagation of waves near any, even fictitious, surface, and it is 
valid also near the interface. With it we do not speak about incident and reflected waves. 
We look for all possible solutions near the interface and select those which correspond to 
our physics. Solution of Eq. (14), gave us wave modes and their speeds, solution of (78) 
will give wave modes and their q or k±. We select in between them, say, one wave with 
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positive and two waves with negative q. They correspond to the incident and reflected 
waves. And we find such a superposition of these waves that satisfies the boundary 
condition (74). Thus we obtain the result absolutely equivalent to that obtained with Eq. 
(14). In the case of a real k± the value of q is tgy?, but q can be also defined for arbitrary 
complex k±, and this is the benefit of the Eq. (78) 

All the waves, incident, reflected or surface ones, should satisfy this equation for the 
given value of T, which is determined by the grazing angle tp of the incident wave and by 
direction of the vector a. Polarization vector A is represented as 

A = an + (3t. (79) 

To find a and [3 we multiply both parts of Eq. (78) by n and r, and obtain a system of 
two linear homogeneous equations, which has solution, when its determinant is equal to 
zero. The resulting equation is a polynomial of the 4-th order in powers of q, and it has 
4 roots. 

For instance, calculations for cosy? = 0.3, which is below 0.5, give all the roots q to 
be real. Two of them are positive: q-y = 3.18, g 4 = 2.47; and the other two are negative: 
q2 = —3.6, qs = —1.4 (such numerations of the roots is for further convenience). These 
roots determine all the possible waves near the surface. The positive roots correspond 
to waves incident on the surface, and the negative ones correspond to waves going away 
from the surface. The root q\ corresponds to the given cosy? = 0.3 of the incident quasi 
transverse wave A^ 2 \ The root q^ shows that, if the incident wave were quasi longitudinal 
one, its grazing angle would be cosy? = 0.37. The negative roots are related to the 

(2) (3) 

reflected waves: q 2 to the quasi transverse, A R , and q 3 — to the quasi longitudinal, A R ', 
ones. 

When cosy? = 0.6 > 0.5 the two roots, q 3 and g 4 , related to quasi longitudinal waves 
become complex conjugate: g 3i4 = 0.32 =f 0.77i We have to take into consideration 
only g 3 , which is related to quasi longitudinal surface wave. This root has correct sign 
of the imaginary part, which warrants an exponential decay of the wave away from the 
interface in the half space z < 0. However it does contain also a real part, which seems to 
make this root unacceptable. In particle physics the wave function ip oc exp(iq'z + q"z) 
at z < means that there is a flux of particles j oc q'exp(2q"z) toward the surface, 
which increases exponentially from z = -co, and shows that during propagation from 
z = — oo toward z = the particles are created from nothing. Intuitively we expect 
the same of elastic waves. However elastic waves behave differently, and because of that 



23 



we call their properties counter intuitive. The numerical calculations of the energy flux 
according to (45) show that the surface quasi longitudinal wave does not create energy 
flux notwithstanding that its k^ has a complex value. 

The most terrible situation seems to occur after cosy? = 0.92. (It is not a fundamental 
constant. It depends on direction of the anisotropy vector a and on values of parameters £ 

(2) 

and E). At some critical y? c i ~ 0.921352 the value of q 2 and therefore of k ± become zero, 
i.e. our anxieties came true! The incident plane wave turns into a surface one! However 
above this critical point the quasi transverse mode does not become of a surface type. Its 
q 2 and therefore k^ do not acquire a negative imaginary part. Instead k^ remains real 
but changes its sign! 

Intuitively we can expect that after reverse of the sign of q 2 the wave becomes propa- 
gating toward the surface. Such a wave should carry the energy also toward the surface. 
Nothing like that! We found that the energy flux of this mode did not change its sign. 
Reflected energy flux related to this mode remains completely equal to the incident flux 
and opposite in direction. It can be understood because the energy flux depends not solely 
on the wave vector k but also on polarization (or oscillation) direction A and anisotropy 
vector a. The direct calculations show that we have no reason to worry nor about energy 
conservation, nor about boundary conditions. They both are satisfied at cosy? > cosy? c i. 

However it is not the end of the story. When we decrease tp below tp c i, the value of 
qi = tgy? decreases and the energy flux of the incident wave decreases too. This is natural. 
Reflected flux decreases in the same way, though q 2 > steadily increases. But there is a 
second critical point y? c2 , where q\ = q 2 , and the energy flux density of the incident wave 
becomes zero! After this point the roles of the two roots q\ and q 2 do exchange. The 
incident wave gives the flux away from the surface, and the reflected wave — toward it. 
Of course it means that the incident wave of the mode does not exist below (p c2 \ All 
that leads us to an interesting conclusion, but before going to it let's discuss the surface 
waves on a free surface in an anisotropic media. 

D. Surface waves 

The first question is: whether the surface waves do exist? From the very beginning it 
was found that if we require that a surface wave to decay away from the surface with a 
real exponent, the equation for the speed of the surface wave leads to a complex value of 
cr, which means that the surface waves are leaky, and therefore cannot be accepted as a 
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stationary solution of the wave equation. However an experience with quasi longitudinal 
surface waves had shown that we can accept a complex exponent. Then we may expect 
to find a real value for cr. 

A surface wave (the Rayleigh one) satisfies the same equation (78) as any other wave, 
but with T = (cr/q) 2 , which we earlier (see Eq. (58)) denoted as x 2 . With it we rewrite 
Eq. (78) as 



The easiest way is to find x by try and error method. We suggest some value of x — x\ < 1, 
seek the solution of (80) in the form (79). Multiply both parts by n and r, obtain two 
homogeneous linear equations for a and f3, find its determinant, which is a polynomial of 
4-th order in powers of q: D±(cj), and find its roots qi (i=l-4). If all the roots are complex, 
we choose two of them with negative imaginary parts, say q 2 and q 3 . For them we find 
a (2 ' 3) , /3 (2 ' 3) and A (2 ' 3) (g 2 ,3)- After that we use (30) and obtain _B (2,3) . With these vectors 
we construct a linear combination, which satisfies boundary conditions 



Multiplication of this equation by n and r gives us again a system of two equations. It 
is resolvable, if its determinant D(xi, q 2 , (73) is equal to zero. 

For an arbitrary chosen x\ the determinant D(x\, 52, 53) 7^ 0. Instead it is a complex 
number, say D(x 1 , q 2 , 93) — yi + iz\. Then we try another x 2 till we find D(x 2 , q 2 , 53) = 
—y 2 — iz 2 , where signs of yi >2 and respectively of zi >2 are the same. After that by narrowing 
the interval x±,x 2 we find the limiting point x , where D(x ,q 2 ,q 3 ) = 0. The Rayleigh 
speed is cr = x c t . In the case of E = 1.5, £ = 0.4 and cos# a = 0.4 we got cr = 0.6066q. 



We formulated the theory of elastic waves in isotropic media with the help of complex 
vector wave functions like in particle physics. We considered reflection and refraction of 
waves at an interface with mode conversion or in other words with double splitting of the 
reflected and refracted waves. We had shown that in the case of a transverse incident wave 
there is a critical grazing angle ip c , below which the longitudinal reflected wave becomes 
of the surface type with a speed in the interval (c t , q). The speed of the Rayleigh wave is 
a root of the equation 1/r = where r is one of reflection amplitudes. 



x 2 - 1 - q 2 + f (fc • a) 2 ] A = Ek(k ■ A) — 
i (o[(l + q 2 ){a -A) + {k- a)(k ■ A)} + k(k • a){a • A 



)])■ 



(80) 





V. CONCLUSION 
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The theory for isotropic media was generalized to anisotropic ones with a single vector 
of anisotropy and a specific term in the free energy of deformation. In such media the 
transverse and longitudinal waves become hybridization, reflection and refraction at an 
interface is accompanied in general by triple splitting of reflected and refracted waves, 
and all the reflected waves are nonspecular. 

In some cases, when speeds of all the reflected waves are higher than that of the 
incident one, a plane wave at some critical grazing angle ip cl can be expected to completely 
transform into a surface one, which violates the energy conservation law. Because of some 
counter intuitive properties of elastic waves in anisotropic media such a transformation 
does not take place. However there are two critical points in the grazing angle of the 
incident wave, which can be considered as a hint that a nonlinearity should come into 
play near these points. If the nonlinearity is included, then the phenomenon, like that one 
shown in fig. 5b), could be possible. Transformation of a plane wave into a surface one 
should lead to an exponential grows of the surface wave amplitude, which can be related 
to such natural phenomena as the devastating earth quakes. 

In many other aspects the wave theory for anisotropic media is a alike to those for 
isotropic ones. It predicts the Rayleigh surface wave on a free surface and the Stonley 
wave on an interface. Theses surface waves have complex normal components of the wave 
vector, however it does not lead to violation of energy conservation, because the real part 
of this normal component does not create an energy flux from the surface. We would like 
to stress that the surface waves, which exponentially decay away from the interface and 
at the same time oscillate, are not so called "leaky surface waves", because their energy 
leaks nowhere. The leaky surface waves cannot exist as a stationary solution of the 
wave equation without introduction of some losses because of nonlinearity or scattering, 
otherwise they violate the law of energy conservation. 
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